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Resumen 


Se deduce una solución analítica de la ecuación de Richards utilizando las hipótesis 
de Green y Ampt para describir la infiltración del agua en un suelo con un manto 
freático somero. Se supone un perfil de humedad inicial lineal, de modo que el 
valor mínimo se encuentra en la superficie del suelo y el máximo en la superficie del 
manto freático. Se acepta una variación lineal de la succión en el frente de saturación, 
de modo que es máxima en la superficie del suelo y nula en el manto freático. La 
ecuación de Green y Ampt se deduce cuando la profundidad del manto freático 
tiende a infinito en la solución. La solución se compara con una solución numérica 
para la condición de contenido de humedad inicial lineal, supuesta para deducir la 
primera, y para una distribución hidrostática inicial de las presiones, con buenos 
resultados en ambos casos. En el modelo de Lewis y Milne se introduce la solución 
de la infiltración para describir tres pruebas de avance en el riego por melgas del 
cultivo de arroz en La Chontalpa, Tabasco, México, a diferentes profundidades 
del manto freático y caudales de riego. En la primera se calibran los parámetros 
relativos a la infiltración y a la ley de resistencia al flujo del agua en la superficie 
del suelo, y en las otras dos se lleva a cabo la predicción de la evolución del frente 
de avance; en las tres pruebas, las curvas teóricas de avance son muy cercanas a 
las curvas experimentales. La solución de la infiltración establecida se utiliza en el 
diseño del riego por gravedad en suelos con un manto freático somero, como ocurre 
en zonas agrícolas con riego o con drenaje lento. 


Palabras clave: perfil de humedad inicial lineal, variación lineal de la succión en el 
frente de saturación. 


con soluciones numéricas o analíticas de la 
ecuación de Richards (1931). 
Entre las soluciones analíticas de la infiltra- 


Introducción 


En el estudio del movimiento del agua en el 


riego por gravedad se han utilizado ecuaciones 
racionales para describir el movimiento en 
la superficie del suelo y la infiltración. El 
movimiento en la superficie se modela gene- 
ralmente con las ecuaciones de Saint-Venant 
(1871) y sus simplificaciones, como los modelos 
de onda difusiva, de inercia nula e hidrológico 
(Woolhiser, 1975; Strelkoff y Katopodes, 1977; 
Lewis y Milne, 1938). La infiltración se modela 


ción se encuentran la de Philip (1957, 1969), 
Parlange et al. (1982) y Fuentes et al. (2001), entre 
otras. Estas soluciones se obtienen suponiendo 
una columna semi-infinita de suelo con un 
contenido de humedad inicial constante a lo 
largo de la misma. Esta condición no es aplicable 
en suelos con mantos freáticos someros, ya que el 
perfil de humedad no es uniforme, como ocurre 
en zonas agrícolas de riego. 
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El objetivo de este trabajo es la deducción 
de una solución analítica de la ecuación de 
Richards basada en las hipótesis de Green 
y Ampt que incluya la presencia del manto 
freático en el perfil del suelo. 

La solución obtenida se compara con una 
solución numérica de la ecuación de Richards 
y se utilizará en la descripción de la fase de 
avance en el riego por gravedad en el cultivo 
de arroz en La Chontalpa, Tabasco, México, a 
través del modelo hidrológico propuesto por 
Lewis y Milne. 


Teoría 
Ecuaciones de base 


La ecuación de Richards se utiliza ampliamente 
para estudiar la transferencia de agua en 
el suelo. En su forma unidimensional en la 
transferencia vertical se escribe como sigue: 


dy _0 dy | dKoy 
car <A ULER 
SE | 06] a Y 


la cual resulta de la ecuación de continuidad 
deducida a partir de la aplicación del principio 
de conservación de la masa: 


%__% (2) 


y de la ley de Darcy generalizada a los suelos 
no saturados: 


donde t es el tiempo; H, el potencial hidráulico 
igual a la suma del potencial de presión y y 
del potencial gravitacional z asimilado a la 
coordenada vertical orientada positivamente 
hacia abajo; q, el flujo; K(yw), la conductividad 
hidráulica como una función de la presión; 0, 


el contenido volumétrico de humedad; 0(w), la 
curva de retención de humedad, y su derivada 
C(w) = 40 /dy es la capacidad específica. 

Se han establecido algunas soluciones analí- 
ticas del fenómeno de la infiltración vertical 
basadas en la ecuación (1), correspondientes 
a formas particulares de las características 
hidrodinámicas y de las condiciones iniciales 
y de frontera. Destaca particularmente la 
solución de Parlange et al. (1982), quienes 
suponen una relación entre la difusividad 
hidráulica definida por D(0) = K(0)dy/d8 y 
la derivada de la conductividad hidráulica 
con respecto al contenido de humedad dK/ 
d8, para una columna semi-infinita de suelo 
con condición de humedad inicial constante 
(0,) y presión nula en la superficie del suelo; 
la solución incluye como casos extremos las 
ecuaciones de Green y Ampt (1911), y de Talsma 
y Parlange (1972). La primera corresponde 
al caso en que D(0) tiene el comportamiento 
de una densidad de Dirac y K(0) es continua, 
humedad, 
mientras que la segunda cuando tanto como 


proporcional al contenido de 
D(0) y dK/d8 presentan un comportamiento 
proporcional a dicha densidad. Parlange et al. 
(1982) muestran que su solución es aplicable 
para describir la infiltración en suelos de 
textura contrastantes, desde una arena hasta 
una arcilla. El efecto de una carga positiva lo 
introducen Parlange et al. (1985), aceptando 
la hipótesis de Green y Ampt en cuanto a la 
presencia de una zona totalmente saturada en 
el perfil de humedad; una zona saturada en 
tensión la agregan Haverkamp et al. (1990). 

En otra línea de investigación, Fuentes et 
al. (2001) ofrecen una solución que, para las 
mismas condiciones iniciales y de frontera, 
suponen que D(8) es constante y que dK/d8 es 
lineal, generalizando el término convectivo de 
la ecuación de Burgers (1948) planteada en un 
problema de aerodinámica y de la cual se han 
construido soluciones analíticas. El resultado 
incluye la solución lineal clásica y la presentada 
por Knigth (1973). 

Las soluciones mencionadas no se pueden 
aplicar en suelos con perfil de humedad inicial 
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no uniforme, sin cometer posiblemente errores 
importantes, como el que se presenta cuando 
existe un manto freático somero en el perfil 
del suelo. Se va a mostrar que la presencia 
del manto puede incorporarse de una manera 
directa en la ecuación de Green y Ampt, razón 
por la cual se expone la solución original. 


La infiltración según Green y Ampt 


La ecuación de Green y Ampt se establece con 
las siguientes hipótesis: 


a) El contenido de humedad inicial (0) es 
constante a lo largo de la columna de suelo, 
es decir8,= 0, enO < z. 


b) Durante la 
zonas de humedecimiento, una totalmente 


infiltración se forman dos 


saturada, es decir8=09en0 < z < z y donde 
9, es el contenido de humedad a saturación 
y z/t) es la posición del frente saturado, 
y Otra seca con el contenido de humedad 
inicial, es decir O = 6, en 2,< 2. La primera 
zona evoca un pistón (figura 1). 
c) La distribución de las presiones es 
hidrostática en la zona saturada, siendo y 
Más 


impuesta en la superficie del suelo, y y = y, 


en z=0, donde hy 2 0 es la presión 


en z = 2, donde y,= -—h,< 0 es una presión 
negativa efectiva en el frente de saturación 
y O <h, es la succión correspondiente. Lo 
anterior se traduce en y = E Mo, + h)z/ 
Éd 

La segunda hipótesis implica que 99/0t = 0 
en0O<z< 2,Y de acuerdo con la ecuación (2) se 
tiene 0q/0z = 0, esto quiere decir que el flujo 
de Darcy es independiente de z en la zona 
saturada y es igual al flujo de infiltración que 
ocurre en la superficie qs(t) = dl/dt, donde 
I(t) = (0.-0,)2/t), de acuerdo con la primera 
hipótesis, es la lámina de agua infiltrada 
acumulada. La tercera hipótesis proporciona el 
gradiente hidráulico 0H /dz =-[1 + Mo, +h ) /z al 
De la ley de Darcy, ecuación (3), se deduce: 


Oo Oi Os 


V 


Figura 1. La infiltración en un suelo con manto freático y 
perfil de humedad inicial lineal, aceptando las hipótesis de 
Green y Ampt (1911). 


dl _ Bop HN . e 
2-5, 70 ; 1(t) =A8z, (+) (4) 


donde A9=0 —6, y K.=K(0) es la conductividad 
hidráulica a saturación. 

La integración de la ecuación (4), sujeta a la 
condición 1] = 0 en t =0, conduce a la ecuación 
de infiltración de Green y Ampt, para h,,, 
independiente del tiempo: 


I 


En una columna de suelo con un manto 
freático localizado a una profundidad P,, la 
ecuación (5) podría aplicarse hasta el tiempo t, 
que corresponde al momento en que el frente 
de saturación alcanza el manto freático, con una 
lámina infiltrada acumulada máxima ] E A0P,. 
Sin embargo, su aplicabilidad es cuestionable, 
ya que la presencia de un manto freático induce 
un perfil de humedad no uniforme (Gardner, 
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1958; Kunze et al., 1985) y la condición inicial 
depende de la profundidad 0 (2), de modo que 
9,(0) = 6, y 0/(P) =0 


La infiltración en una columna de suelo con manto 
freático 


Se denota por 0(z) la distribución inicial 
de la humedad; su desarrollo en potencias 
alrededor de z = 0 es 0 (2) = 0, + 0,(0)z + O(23. 
Si se retienen los términos hasta el lineal, 0,(0) 
debe ser reemplazada por (0,- 0,)/P, con la 
finalidad de satisfacer la condición 0(P) = 0. 
Así, la primera hipótesis de Green y Ampt se 
reemplaza aquí por una dependencia lineal del 
contenido de humedad inicial con respecto a la 
profundidad: 


9; (2) Ue (0,- 0) (6) 
f 


la cual contiene el caso de un contenido de 
humedad constante haciendo P, > co. En la 
figura 1 se muestra un esquema de la infiltra- 
ción según Green y Ampt, con el contenido de 
humedad inicial definido por la ecuación (6). 

La lámina infiltrada acumulada en un flujo 
en pistón con contenido de humedad inicial 
variable está definida por: 


zy(t) 


I()= | [0,-0,(2)]az (7) 


La introducción de la ecuación (6) en la 
ecuación (7) proporciona la relación siguiente 
entre la lámina infiltrada y la posición del frente 
de saturación: 


1(0/M=1-[1-2,(0/2,]: 


donde !,, es la lámina infiltrada máxima. 


In=340P, (8) 


La succión en el frente de saturación 
depende de los contenidos de humedad inicial 
y saturación, es decir h, = h0,0.). Cuando el 
frente de saturación está en la superficie se 


tiene el valor h (0 /0.), y cuando el frente de 
saturación está en la superficie del manto 
freático el valor de la succión es nula; su 
desarrollo en potencias alrededor de 0,= 0, es 
h/(0,8,) = h(0,,0,) + (H,/98)o;.., (0,= 0,) + O[0,- 
9 y]. Si se retienen los términos hasta el lineal, 
(9h,/ 00), se debe reemplazar por -h(0,,0,)/ 
(0. — 0,) para satisfacer la condición en el 
manto freático: h,(0,0,) = 8,831 - (0,- 0,)/ 
(0.8 )]. Utilizando la ecuación (6) se obtiene la 
dependencia lineal siguiente: 


ny[0,(2p), 0,] = 11, (0,, 0)[1- a (9) 
a 


La ecuación (4) se transforma en la ecuación 
diferencial siguiente: 


All, y rote tti120|, 


ej 
2=P,[1- 1] 


donde se ha escrito h,en lugar de (0,8), para 
que corresponda al parámetro de la ecuación 
original de Green y Ampt. Se debe notar que 
la velocidad con la que llega el frente al manto 
freático (2,= P) es q,=K(1 + Mo. l PO, en lugar 


de q,= K; esta última situación ocurre cuando 


La integración de la ecuación (10) con la 
condición 1 = 0 en t = 0 conduce a la ecuación 
de infiltración que sigue, cuando h,+* P;: 


P; T- 2Py (.,, +24) SID 


(11) 


¿22 ta ly) Ls (h 


Cn la) 


y cuando h,= P,se tiene: 
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Se debe notar que la ecuación (11) se 
reduce a la ecuación de infiltración de Green y 
Ampt, ecuación (5), cuando P, > «o. En efecto, 
considerando que ] ¿E A0P, en este límite se 
tiene 1 — (1 - 1/1,)P= 1/21 ¿=1/48P, El tercer 
término es del orden de 1/P ¿Y tiende a cero; en 
el segundo término, el argumento del logaritmo 
tiende a 1 + I/A0(h,., 
A0(h,,,, +Hh); 


sup 


+ h) y su coeficiente a 
finalmente, en el primer término, el 
coeficiente de I tiende a la unidad. Tomando en 
cuenta la definición A = A0(..., + hp, se deduce 
la ecuación de Green y Ampt. 

En tiempos cortos, la solución se desarrolla 
en la raíz cuadrada del tiempo: 


I=S Vt+At+0 (4%) (13) 


donde: 


ge _2 Moup + 3h 
Si =2K, (M1; +11.) A8; A K, (14) 


El parámetro S, se denomina la sorbilidad 
total y se infiere que su cuadrado es el resultado 
de las propiedades del suelo y de la condición 
en la frontera; es decir, S? = NE donde 
5”=2K h,A8 y S”,,=2K h,,, AO. El parámetro 
S es propiamente la sorbilidad del suelo (Philip, 
1957). También cuando P,> «o se obtiene el 
resultado clásico A= 3K, del desarrollo de la 
ecuación de Green y Ampt. 

Para deducir el comportamiento cuando 
t>t,y dondet,,eseltiempo correspondiente 
a la lámina infiltrada máxima, primero se 
deduce de las ecuaciones (11) y (12): 


P; L a (1,,¿+P,) (a, + La 
3 
có: (E, =14) 
h +P:K,t El 
= 2P4 (1... + hy 
m0 e] 2 
sat ly (b- ll 


P¿l, 
h¿=P,: K.tr= M (16) 
f f sM 
3 (Maup+ +) 
Restando las ecuaciones (11) y (15), 


y las ecuaciones (12) y (16), se deduce el 
comportamiento buscado: 


hpePy: 1-1y+K,(1-hs/P, ME tm) (17) 


hy=Py: I-Iy+K, (L+Mop/Pp)(- tm) (18) 


Estimación de parámetros 


El parámetro h, se calcula a partir de las 
características hidrodinámicas del suelo, utili- 
zando la expresión que define la sorbilidad: 


¡0(w)-0 


£[9()] 


"K(wdy (19) 


ES -aj* oO e. 


donde f(0) es la relación de concentración 
de flujos y 0, = 0(y.,), siendo y, la presión 
correspondiente a 0. 

La expresión que resulta para la succión en 
el frente de saturación es: 


| 20 ay; 09-202 0) 


Few] 9,-0, 


donde K(w) = K(w)/K, es la conductividad 
hidráulica relativa a su valor en saturación. 

La relación de concentración de flujos para 
una difusividad de Dirac es (O) = O. Mediante 
un procedimiento de optimización, Parlange 
(1975) recomienda: 


1)-= 


21 
o (en 

El contenido de humedad inicial en la 
superficie se estima a partir de la curva de 
retención de humedad. En efecto, Saucedo 
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et al. (2003) muestran que el supuesto de una 
distribución hidrostática de las presiones en el 
tiempo inicial rinde cuenta de la evolución del 
frente de onda en la fase de avance en el riego 
por gravedad en suelos con mantos freáticos 
someros. Este supuesto permite reemplazar el 
contenido de humedad inicial en la ecuación 
(7) por la curva de retención 9(w), considerando 
que y = z — P,. Para la lámina de infiltración 
máxima se tiene: 


P. 
lu= f[o.- 9(2-P,)]az (22) 
0 


Así, la ecuación (22) permite estimar un 
contenido de humedad inicial en la superficie, 
de modo que se conserve la masa total. De 
acuerdo con la ecuación (8), LS A0P,, se 
tiene: 


1 
Ps 


Pe 
2A0 = fo. 0(2 —Py)] dz (23) 
0 


si 0,< 8, entonces 0 ,=8,, donde 6, es el contenido 
de humedad residual de Brooks y Corey 
(1964), definido de modo que la conductividad 
hidráulica asociada es cero. 


Comparación entre las soluciones analítica y 
numérica de la infiltración 


Para comparar la solución analítica con la 
solución numérica de la ecuación de Richards 
se utiliza un suelo arcilloso de La Chontalpa, 
caracterizado hidrodinámicamente por Sau- 
cedo et al. (2003). Estos autores representan 
la curva de retención con la ecuación de Van 
Genuchten (1980) y la curva de conductividad 
hidráulica con la ecuación de Brooks y Corey 
(1964), como es sugerido por Fuentes et al. 
(1992): 


9 (w)=0,+(8,-0,) mL ; m=1-2/n (24) 


K(0) =K, E (25) 


9,0, 


Los valores de los parámetros en este suelo 
arcilloso son 0,= 0.5245 cm*/cm', O, =0 cm*/cm*, 
y,=-15 cm, K =1.86 cm/h, m= 0.066, n =21.19. 

Se utiliza el código numérico propuesto 
por Zataráin et al. (1998) para obtener 
la solución de la ecuación de Richards, 
ecuación (1), sujeta a la condición inicial 
definida por la ecuación (6) con 9,= 0.4 cm*/ 
cm? y diferentes profundidades del manto 
freático P,= 50, 100, 150 cm. La presión en 
la superficie seleccionada es h,,, = 10 cm, 
que es del orden de magnitud del tirante de 
agua formado en el riego por melgas en el 
cultivo de arroz. En la figura 2 se muestra la 
evolución de la lámina infiltrada calculada 
con la ecuación de Richards y con la solución 
analítica propuesta; el buen acuerdo entre 
ambas soluciones es evidente. 

En la figura 3 se tiene la misma compara- 
ción que en la figura 2; pero en este caso, la 
ecuación de Richards se resuelve suponiendo 
una distribución hidrostática de las presiones. 
El acuerdo sigue siendo bueno, aunque la 
solución analítica se obtiene para una condición 
inicial diferente. 


— Richards 
=—- Solución analítica 


PF 50 cm 


Lámina infiltrada (cm) 


0 Y Y T Y T 
0 50 100 150 200 250 300 
Tiempo (min) 


Figura 2. Comparación entre la solución numérica de la 
ecuación de Richards sujeta a un perfil de humedad inicial 
lineal y la solución analítica definida por la ecuación (11). 
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E — Richards 
12 noe das 
El -—- Solución analítica 
"o > 
S 1073 
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E 6 Pf =100 cm 
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A 

2 +5 

0 s 

0 100 200 300 400 
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Figura 3. Comparación entre la solución numérica de la 
ecuación de Richards sujeta a una distribución hidrostática 
inicial de las presiones y la solución analítica definida por 
la ecuación (11). 


Aplicaciones 


Una de las posibles aplicaciones de la ecuación 
de infiltración deducida es la descripción de 
la evolución del frente de onda en la fase de 
avance del riego por gravedad en suelos con 
mantos freáticos someros. 

Las pruebas de riego por melgas del cultivo 
de arroz fueron realizadas por Pacheco (1994) 
y modeladas mediante el acoplamiento de las 
ecuaciones de Saint-Venant y de Richards por 
Saucedo et al. (2003). Por simplicidad, aquí se 
utiliza el modelo hidrológico de Lewis y Milne 
(1938): 


xp(t) 


qot=hxy(0)+ | I(t-+,)dx (26) 


donde q, es el caudal de aporte de agua por 
unidad de ancho de la melga o caudal unitario 
de aporte; h, un tirante de agua medio en la 
superficie del suelo; 1 = t — £, el tiempo de 
contacto del agua en un punto x de la melga, 
y es igual a la diferencia entre el tiempo de 
avance (t) y el tiempo en que la onda arriba al 
punto en cuestión (t). 

La ecuación (26) puede ser resuelta 
numéricamente. La solución para At constante, 
utilizando la regla de los trapecios, es la 


siguiente: 


qot¡—Hx;x 5 Lat Mx, 
X= Xj4 + n +0 en 
1=2,3, 
con t,=At, t,=jAt, x,=0 y x,=,t,/(h + c,1,). El 


coeficiente de forma está en L< c;<1; se puede 
tomar el valor c, =2. 

Para el cálculo del tirante medio en función 
del caudal unitario de aporte es necesario 
introducir una ley de resistencia al flujo. Se 
utiliza la ley de resistencia en potencia siguiente 


(Fuentes et al., 2004): 


d 
y? 
ES 
vV 


donde q es el caudal unitario; h, el tirante 
de agua; J, la pendiente de fricción; v, el 
coeficiente de viscosidad cinemática; g, la 


(28) 


aceleración gravitacional; Lb<d<l, el límite 
inferior corresponde a la ley de Chézy (1776) 
y el superior a la ley de Poiseuille (1841); y k 
es un coeficiente adimensional. Se suponen los 
valores v = 1x10* m?/s y g = 9.81 m/s?. 

A partir de un análisis de la forma de la 
onda de avance en tiempos largos, Zataráin et 
al. (2003) concluyen que el tirante medio del 
modelo hidrológico puede ser estimado con la 
fórmula: 


138dy 2 y13 
mt dl (6 8) 
5d+1 kv 8). 


donde h, es el tirante normal y /, es la pendiente 
topográfica de la melga. El tirante medio 
es tomado como el valor de la presión en la 
= Hi 

Los caudales unitarios de las tres pruebas 


superficie del suelo: /,,. 


de riego son reportados en el cuadro 1; la 
pendiente topográfica de la melga es J,=0.00085 
m/m. Tomando el caudal unitario máximo, el 
criterio del número de Reynolds Re «< q,/v=1 428 
indica que el régimen de flujo es laminar, es decir 
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d = 1. En cuanto al parámetro k, Saucedo et al. 
(2003) lo calibran de modo de reproducir la 
curva de avance en la primera prueba de riego; 
con las ecuaciones de Saint-Venant y Richards 
obtienen k=0.0033; las dos pruebas restantes se 
predicen correctamente utilizando este valor. 

Dado que la evolución del tirante 
proporcionada por las ecuaciones completas 
de Saint-Venant se reemplaza por un tirante 
medio y constante en el modelo hidrológico, 
la calibración del parámetro k en la primera 
prueba de riego es diferente; el valor óptimo 
es k=0.0025. En el cuadro 1 se muestran los 
valores de la infiltración máxima calculada 
con la ecuación (22), el contenido de humedad 
en la superficie calculada con la ecuación 
(23), la succión en el frente de saturación en 
la superficie calculado con las ecuaciones 
(20) y (21), y el tirante medio calculado con la 
ecuación (29). 

En la figura 4 se muestran las tres curvas 
de avance. En la primera prueba de riego, la 
curva teórica es de calibración, mientras que en 
las pruebas segunda y tercera, las curvas son 
de predicción. El buen ajuste entre ellas indica 
que la solución analítica planteada puede des- 
cribir la evolución de la lámina infiltrada en un 
suelo con manto freático en su perfil. 

Los parámetros de la solución analítica 
h, y K, se han calculado a partir de las 
características hidrodinámicas obtenidas a 
su vez mediante la calibración de la curva de 
avance con las ecuaciones de Saint-Venant y 
Richards. Sin embargo, los parámetros pueden 
obtenerse mediante la calibración directa de 
la curva de avance con el modelo hidrológico 
y la solución analítica de la infiltración; el 
parámetro h, como se observa en el cuadro 1, 


100 
o Pf> 152 cm 


so] Pf= 52 cm 


A Pf= 50 cm 


Tiempo (min) 
E 


o 2 40 60 80 100 


Distancia (m) 


120 


Figura 4. Posición del frente de avance en el riego por 
melgas en el cultivo de arroz en La Chontalpa, Tabasco, 
México. En el cuadro 1 se reportan los valores de los 
parámetros del modelo teórico. 


es relativamente constante en un buen rango 
del contenido de humedad en la superficie del 
suelo y puede tomarse como el que resulta de 
la primera prueba de riego. Los parámetros se 
estiman también a partir de pruebas clásicas 
de infiltración. 


Conclusiones 


Se ha deducido una solución analítica de la 
ecuación de Richards de la transferencia de 
agua en el suelo utilizando las hipótesis de 
Green y Ampt para describir la infiltración 
del agua en un suelo con un manto freático en 
su perfil. Por una parte, se supone un perfil 
de humedad inicial lineal, de modo que el 
valor mínimo se encuentra en la superficie del 
suelo y el máximo en la superficie del manto 
freático. Por otra parte, se acepta una variación 
lineal de la succión en el frente de saturación, 
de modo que es máxima en la superficie del 


Cuadro 1. Valores de los parámetros en las tres pruebas de riego en el cultivo de arroz realizadas 


en una melga de 10.5 m de ancho y 100 m de largo, en La Chontalpa, Tabasco, México. 


Prueba q, (m*/s/m) P,(cm) 1, (cm) 0, (cm*/cm”) h,(cm) h (cm) 
ll 0.001428 152 14.54 0.3331 14.79 2.73 
2 0.001428 50 2 IS 0.4386 18:55 2.73 
3 0.001238 52 2.32 0.4353 13.64 2.60 
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suelo y nula en el manto freático. La solución 
analítica tiende a la ecuación de Green y Ampt 
conforme la profundidad del manto freático se 
incrementa. La solución se compara con una 
solución numérica de la ecuación de Richards 
para la condición de contenido de humedad 
inicial lineal, supuesta para deducir la primera, 
y para una distribución hidrostática inicial de 
las presiones, con buenos resultados en ambos 
casos. 

En el modelo de Lewis y Milne de la fase de 
avance del riego por gravedad, la solución de 
la infiltración se introduce para la descripción 
de tres pruebas de riego por melgas realizadas 
en el cultivo de arrozen La Chontalpa, Tabasco, 
México, a diferentes profundidades del manto 
freático. En la primera prueba se calibraron 
los parámetros relativos a la infiltración y 
a la ley de resistencia al flujo del agua en la 
superficie del suelo; en las pruebas segunda 
y tercera se llevó a cabo la predicción de la 
evolución del frente de avance. Tanto en la 
calibración como en la predicción, las curvas 
teóricas de avance son muy cercanas a las 
curvas experimentales. 

A partir de lo anterior, la solución de 
infiltración propuesta puede ser utilizada en 
el diseño del riego por gravedad en suelos 
con un manto freático somero, como ocurre en 
zonas agrícolas de riego o suelos con drenaje 
natural lento que inducen tales condiciones. 
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Abstract 


FUENTES, C., CHÁVEZ, C. € ZATARÁIN, E. An analytical solution for infiltration in 
soils with a shallow water table: application to gravity irrigation. Water Technology and 
Sciences, formerly Hydraulic engineering in Mexico (in Spanish). Vol. I, no. 3, July- 
September, 2010, pp. 39-49. 


An analytical solution of the Richards equation is deduced using Green and Ampt hypotheses 
to describe the infiltration in soils with a shallow water table. A linear initial moisture profile 
is assumed, so that the minimum value is at the soil surface and the maximum value at the 
water table level. A linear variation of the driver pressure is accepted, it is maximum at the 
soil surface and zero at the water table level. The Green and Ampt equation is deduced from 
the solution when the water table depth tends to infinite. The solution is compared with a 
numerical solution for linear initial moisture profile and for an initial hydrostatic pressure 
distribution, with good results in both cases. In the Lewis and Milne model, the infiltration 
solution is introduced to describe three tests of advance phase in border irrigation of the rice 
culture at La Chontalpa, Tabasco, Mexico, to different water table depths. In the first one 
the parameters regarding the infiltration and the flow resistance law at the soil surface are 
calibrated and in the other two tests, the prediction of the advance front evolution is carried 
out; in the three tests, the advance theoretical curves are very near to the experimental curves. 
The established infiltration solution can be used designing border irrigation in soils with 
shallow water table like agricultural zones under irrigation or with slow drainage. 


Keywords: linear initial moisture profile, linear variation of the driver pressure. 
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